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We study the equilibration dynamics of closed finite quantum systems and address the question of the time 
needed for the system to equilibrate. In particular we focus on the scaling of the equilibration time Teq with 
the system size L. For clean systems we give general arguments predicting Teq — O (L") for clustering initial 
states, while for small quenches around a critical point we find Teq = O (L'') where ^ is the dynamical critical 
exponent. We then analyze noisy systems where exponentially large time scales are known to exist. Specifically 
we consider the tight-binding model with diagonal impurities and give numerical evidence that in this case 
Teq ~ Be^^ where B, C, ij) are observable dependent constants. Finally, we consider another noisy system 
whose evolution dynamics is randomly sampled from a circular unitary ensemble. Here, we are able to prove 
analytically that Teq = O (1), thus showing that noise alone is not sufficient for slow equilibration dynamics. 
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I. INTRODUCTION 

Consider a quantum system initialized in a given state and 
then allowed to evolve undisturbed under the action of a time 
independent Hamiltonian. Experimentally accessible quanti- 
ties are expectation values of physical observables A at a given 
time {A{t)). The timescale at which such an expectation value 
relaxes to equilibrium identifies the equilibration time of the 
particular dynamics. In infinite systems, equilibration times 
are typically extracted from the exponential decay of observ- 
ables or correlation functions towards their equilibrium val- 
ues. For finite systems, or more generally for systems with 
discrete energy spectrum, the dynamics is quasi-periodic and 
such exponential decay cannot occur, thus necessitating a dif- 
ferent definition of equilibration time. In these cases a mean- 
ingful definition of equilibration time is the first time instant 
for which the value of an observable equals its equilibrium 
value, i.e. given an observable A and its equilibrium value A, 
Toq is the smallest t for which A{t) ^ A\l\. 

In this paper we study the behavior of equilibration times 
according to this definition for various physical systems. In 
particular we are interested in the scaling behavior of the equi- 
libration time as a function of the linear size of the system 

ma. 

We first consider equilibration times in clean systems. Us- 
ing the Loschmidt echo as a particular observable, we are able 
to give general estimates for the scaling of the Teq as a func- 
tion of length. Given a local Hamiltonian H — ^ - hi and a 
sufficiently clustering initial state |'i/'o)^ i-S- gi.j = {hi^j) ^ 
{hi) (hj) decays sufficiently fast as a function of i — j at large 
separations, we find that Toq — 0{L^), i.e. the equilibration 
time is independent of the system size. This scaling becomes 
Teq = 0{L'^) where C is the dynamical critical exponent in 
case both the initial state and the evolution Hamiltonian are 
close to a critical point. These general findings are checked 
explicitly for the Ising chain in a transverse field. However, 
one should be cautious that these arguments may need mod- 
ifications for other observables e.g. such as those undergoing 
spontaneous symmetry breaking. 



We then turn our attention to systems with random impuri- 
ties. According to intuition based on classical models, one 
expects, in general, a slower transient approach to equilib- 
rium in noisy systems compared to clean systems. Indeed, 
it is known that exponentially large time-scales are present 
in glassy systems ||2l. To check and further understand this 
conjecture, we compute numerically the equilibration time 
for the tight-binding model with diagonal impurities, some- 
times called Anderson model ([31. A similar model (albeit with 
pseudo-random diagonal elements) has recently been studied 
in |4|, where it was found that observables relax following 
a power-law behavior. Such power-law equilibration pattern 
(observed also in |5|for another noisy system) is itself a sig- 
nature of slow equilibration. However, according to our defi- 
nition, equilibration times also depend on the long-time equi- 
librium value. Our findings confirm a very slow equilibration 
dynamics characterized by equihbration times diverging ex- 
ponentially with the system size. 

Disorder alone might not be sufficient to guarantee the pres- 
ence of exponentially large relaxation timescales. To illus- 
trate this point, we analyze a second noisy system. In this 
case the evolution operator is related to a unitary matrix sam- 
pled from the circular unitary ensemble (CUE). This model is 
similar to the ones previously considered in ^JT\, where, ac- 
cording to a different definition, equilibration times decreas- 
ing algebraically with the size were predicted |7|. Using our 
definition we prove analytically that, for the Loschmidt echo, 

Teq = 0(1). 

The paper is organized as follows. In Section II, we de- 
scribe a clean system and the observables and the quench con- 
sidered to study the nature of equilibration. Our numerical 
and analytical results for equilibration timescales and the scal- 
ing of these timescales with system length are outlined for the 
cases of the tight-binding model in Section III and the CUE 
based evolution in Section IV, respectively. We present our 
conclusions in Section V. 
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II. EQUILIBRATION IN CLEAN SYSTEMS 

Before considering noisy systems, let us recall some el- 
ementary yet important facts regarding unitary equilibration 
in clean systems. The system is initialized in some state po 
which evolves unitarily via p{t) = e~**^poe**^. First note 
that, because of the unitary nature of the evolution, p (t) does 
not converge in the strong sense as t — > oo. This is true 
irrespective of the Hilbert space dimension, i.e. also in the 
thermodynamic limit. One can then consider the possibil- 
ity of a weaker convergence by looking at "matrix elements" 
A (t) — tr [Ap (t)] where A is an observable. In the thermo- 
dynamic limit the spectrum becomes continuous and one can 
have limit A{oo) — limt_j.oo -4 (t) for some observables A 
(or appropriately rescaled observables) essentially as a conse- 
quence of Riemann-Lebesgue lemma (see e.g. [81 and also f9\ 
for a recent discussion). For finite systems however, A {t) is a 
trigonometric polynomial and hence, once again, does not ad- 
mit an infinite time limit. For the same reason though, the time 

average A := limT--j.oo Jq A (t) dt exists and is finite. Such 
a time average coincides with the infinite time limit when the 
latter exists, i.e. A — A (oo). So the time average can be seen 
as a mathematical trick, reminiscent of Cesaro summation, to 
obtain the infinite time limit when the function oscillates. Al- 
ternatively the time average can mimic the actual measure- 
ment process. In this case T is the "observation time" which 
one may argue to be very large compared to the time scales of 
the unitary dynamics (see e.g. lilOJ ). Clearly one may want to 
investigate the effect of a finite T, here for simplicity we will 
always take T — cx). 

Because of the above considerations, the standard definition 
of equilibration time, extracted from the exponential decay of 
some time dependent observable or correlation function, does 
not make sense for finite systems, because the dynamics is al- 
most periodic and no exponential decay is possible. Instead, 
in finite systems, expectation values A [t] start from a value 
A (0) which retains memory of the initial state po, and then, 
after an equilibration time, approach an average value A and 



start fluctuating around it with fluctuations 5 A = \J A"^ — A 
due to the finite dimensionality of the system. Clearly the pre- 
cise concept of equilibration time is to some extent arbitrary, 
and different definitions are possible. Ours is the following: 
Tcq is the first time for which A (t) equals the average A, 
i.e. is the first solution of A (Tcq) — A (see, for example, 
fig. [T]). This definition is both simple to implement and physi- 
cally clear. 

Some comments are in order: i) Clearly the precise numer- 
ical value of Tcq is irrelevant, whereas the important infor- 
mation is contained in the scaling dependence of Tcq on the 
system size, Tcq (L). ii) In principle A (t) could intersect A at 
a first time Ti , deviate considerably from the average and in- 
tersect A again at T2, and possibly have multiple intersections 
up to T„ before fluctuations of order SA start to set in. In this 
situation it seems that the equilibration process cannot be cap- 
tured by a single Tcq but rather consists of many time scales. 
In all the situations encountered in our analysis, however, we 
found that equilibration could always be captured by a single 
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Figure 1 : Relaxation of nt (t) — E [Ne {t) /L] with time for lengths 
L = 96, at filling u = N/L = 0.25 and a = £/L ^ 0.5. The 
random distribution has variance a = 0.6 (W = 0.6^3 ~ 1.03). 
The ensemble average is computed using 2000 realizations at each 
time. 



time-scale, iii) The above definition depends on the specific 
observable A and different observables may have, in princi- 
ple, different equilibration times. However, we expect that 
for large classes of observables the scaling dependence with 
L will be the same, iv) Based on the experience with infinite 
systems, where observables decay towards their equilibrium 
value, one could define Tcq as that time for which A (t) is off 
by a small amount from its average value A. This is, after all, 
the definition we use in the thermodynamic limit when the ap- 
proach to equilibrium is exponential. A definition of this sort 
has been used for instance in |7|. In this way, however, Tcq 
depends in principle on the way we define small. 

To test the equilibration time we consider a particular ob- 
servable, the Loschmidt echo (LB) which we define as 



^(0 = |(V^o| 



-itH 



(1) 



where ji/io) is the state in which we initialize the system and H 
the evolution Hamiltonian. The Loschmidt echo has been first 
introduced in the context of quantum chaos (see e.g. [ 1 1 1), and 
is generally given a more general form in that context. Equa- 
tion ([T]i can be seen as the time evolved expectation value of 
the observable \ipo){ipo\, and is also known as survival prob- 
ability. Here we consider the LE because it is amenable of 
a cumulant expansion which correctly approximates £ (t) for 
sufficiently large times of the order of Tcq. This conclusion is 
based on numerical experiments on the Ising model in trans- 
verse field m. A hand waving argument is the following: 
since £ is morally a product of L'^ terms (see below) it should 
be clear that Taylor expansion of In C works better than the 
Taylor expansion of C itself. The cumulant expansion of 
Eq. ([T]i reads 



£ {t) = exp 



. n=l ^ ' 



(2) 



where {■)(, stands for connected average with respect to |-0o) 
Truncating Eq. (j2]i up to the first order we obtain C (t) c: 
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exp [-^2 (Aff2)] (A7J2 ^ l^JJ2^^ _ ^ff'^2y Equating the 

short time expansion to the average value C we get the fol- 
lowing expression for the equilibration time; 



-ln£ 



(3) 



As we will see, the above estimate for T^q works well for 
the Loschmidt echo Eq. ([T}. Some comments are in order 
at this point. First of all, the equilibration time in Eq. ([3]) is 
inversely proportional to the square root of an energy fluc- 
tuation. This is not simply the inverse of an energy gap as 
one might guess naively. Secondly, we would like to compare 
Eq. Q, with another estimate of equilibration time which ap- 
peared recently in a single body setting 1721 . The estimate of 
[12 1 reads Tcq ^ l/(A£',nin)ave, where AEj^m is the min- 
imum energy gap averaged over an energy shell around the 
initial energy (H) |29|. So, apart from the order of aver- 
ages, the equilibration time in IJ2J is inversely proportional 
to a standard deviation of an energy fluctuation as much as in 
Eq. ([3]l. The definitions differ in the numerator which takes 
into account the many-body nature of the problem. Thirdly, 
the estimate Eq. Q is valid only for the equilibration time 
of the Loschmidt echo, and in principle, different observables 
might equilibrate with different time scales. 

We will now provide arguments to estimate Eq. ([3]) which 
first appeared in |8|. For a local Hamiltonian H ~ J^i 
and a sufficiently clustering initial state I^jq), i.e. ^ = 
{hihj) — (hi) (hj) decays sufficiently fast as a function of i—j 
at large separations, all the cumulants in Eq. (j2]) are exten- 
sive in the system size, that is in a d-dimensional system of 
linear size L, {H'^")c oc L'^. This means that at leading or- 
der in L, C (t) ~ e-^*^*^-^ ^30\ and so, by Jensen's inequality 
e^^ < C < exp[i'' maxt / (t)], showing that C is expo- 
nentially small in the system volume (note that we must have 
/ (t) < 0). For this reason it is sometimes useful to consider 
the logarithm of the LE {t) = \n£ (t). The equilibration 

time for would then be given by = ^ —In C/ AH"^ . 

Now, since ln£ ~ f {t)L'^, we see that the equilibration 
times for C and T = \n£ are expected to give the same 
scaling with respect to L. From equation (j2]l we get then 
Teq = 0(1): the equilibration time is independent of the 
system size. This argument fails when one considers small 
quenches close to a quantum critical point as in this case the 
clustering properties of the initial state tends to break down. 
In this case |-0o) is the ground state of Hq ~ H (Aq) where 
the external parameter Aq is close to a quantum critical point. 
One then suddenly changes the parameters by a small amount 
'^o ^ -^0 + 5\ and the system is let evolve undisturbed with 
Hamiltonian H = H (Ao + 5\). In the very small quench 
regime, roughly 5\ ^ min{L^^/'', L^^/^j where v is the 
correlation length exponent, perturbation theory is applicable 
(see below) and the average LE reduces to C ~ | (V'olO) |^, 
where |0) is the ground state of H |13 |. The scaling proper- 
ties of the fidelity |('i/'o|0)| close to a quantum critical point 
have been studied in |14| where it was shown that |(?/'o|0)| ~ 
1 — const. X (5A^-L2(''+^^'^), where C, is the dynamical critical 



exponent, and A the scaling dimension of the operator driv- 
ing the transition. Using similar scaling arguments one can 
show that the variance scales as {H^)c ~ L^Crf-A) jg| fj^gsg 
results are valid in the perturbative regime where KV'olO)] is 
not too far from 1, i.e. 5\ ^ min{L^^/'', L^^/dj. where 

1/ = {d + ( — A)^^ |14|. From equation ^ we now get 
Tcq = O {L''), i.e. the equilibration time diverges for large 
systems. Such a divergence is reminiscent of the critical slow- 
ing down observed in quantum Monte Carlo algorithms. 



A. Quantum Ising model 

As a concrete example we now consider equilibration times 
and in particular the prediction Eq. (|3]l for the quantum Ising 
model undergoing a sudden field quench. The Hamiltonian is 



(4) 



and periodic boundary conditions are used. The system is ini- 
tialized in the ground state of Eq. ^ with parameter h^^\ At 
time t — 0, the parameter is suddenly changed to /i^^) ^nd 
the state is let evolve unitarily with Hamiltonian H (/i2)- The 
model in Eq. Q has critical points in the Ising universality 
class at /i = ±1 with d — v — ( = 1, separating an ordered 
phase \ h\ < 1 from a disordered paramagnetic region \h\ > 1. 
For < 1 the order parameter (erf) becomes non-zero, thus 
breaking the Z2 symmetry (erf — —erf) of the Hamiltonian. 

According to Eq. Q and the discussion of the previous sec- 
tion we expect the following behavior for the equilibration 
time as a function of the quench parameters h2, 6h = — hi 
and system size L: 



oc 



const. 



for /i2 = he 
otherwise 



and Sh <^ L 



(5) 



As a first test we check whether Eq. (|5]l is satisfied for the 
Loschmidt echo itself. The LE for a sudden quench has been 
computed in |15 | (superscripts refer to to different values of 
the coupling constants ft,'*)) 



l-sin^ {S'dk)sin' 



k>0 



(Aft/2 



where = 2-^/1 + h^ + 2hcos{k) is the single particle dis- 



(2,) ± 'I 

persion, S^k = ^k' - ^k' and d'}!' are the Bogoliubov an- 
gles at parameters /i'*) , and the quasimomenta are quantized 
according to fc = tt (2n + 1) /L, n — 0, 1, . . . L/2 — 1 (see 
jSJ for further details). In Fig.|2]we plot the equilibration time 
of the Loschmidt echo versus size for different quench 
parameters computed exactly by solving numerically for the 
first solution of L [t) = L. Indeed Eq. (pjl is satisfied to a high 
accuracy. Moreover the transition between the two behaviors 
of Eq. (j5]l appears to be very sharp. From the numerical anal- 
ysis the following behavior for T^^ emerges valid outside the 
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Figure 2: Equilibration times Tgq for the Loschmidt echo as a func- 
tion of the system size L for different quench parameters h\a 
computed by solving numerically for the first solution of L (t) 



Figure 4: Equilibration times for the transverse magnetization 
(t) = (cti (t)) as a function of the size L for different quench pa- 
£ rameters hi.2. is computed exactly solving for the first solution 



Tcq is 



of (t) = 




220 



E 

h5 




200 



Figure 3: Equilibration times for the Loschmidt echo close to the 
crossover region 5h ~ cL~^ . 



crossover region 5h « cL 




for h2 = hf,, and Sh <^ cL ^ 



^jj^ otherwise 



where the constant c in principle depends on hi, h2, but for 
small quenches close to the critical point tends to c ~ 6.1. 
The typical behavior in the crossover region Sh ~ cL^^ is 
depicted in Fig.|3] 

It is natural to ask whether the predictions of Eq. Q are 
satisfied for observables other than the Loschmidt echo. To 
this end, we consider the transverse magnetization (t) = 
(erf (t)) which can be computed as |8, 16J 

(t) = i ^ cos ') cos {Si)k) 

k 

+ sin {d^i^'^^ sin {5dk) cos (tA^^^) (6) 

From eq. (|6]) we extract the equilibration time from the so- 
lution of {t) = = J2k (^r^) (S'&k)- The 



Figure 5: Equilibration times TJq for the transverse magnetization 
close to the crossover region Sh ~ c'L^^ . For this values of param- 
eters c' ~ 13. 



numerical results, shown in Fig. |4] confirm that Eq. ([3]l is ap- 
plicable to the transfer magnetization too. The numerical re- 
sults can be summarized as 



rpm 

cq 



iSh 



for h2 = he 
otherwise 



and Sh <C c'L 



15 for small quenches close 



with constant given now by c' 
to criticality. 

Finally, let us comment on the approach to equilibrium of 
the order parameter (t) — {af (t)). We recall that, since a 
non-zero ni^ breaks the symmetry of the Hamiltonian, 
must be computed as the clustering part of an equal time 



correlation function: 



As 



usual in symmetry broken phases this requires the thermody- 
namic limit to be taken first, but a finite size approximation for 
systems with periodic boundary conditions can be obtained 



by considering the correlation at half chain [m] 

(^f (i)<+L/2 it))- 
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We do not expect formula ^ to reproduce the equilibration 
time correctly for the order parameter, as it does not distin- 
guish whether we are in the ordered phase or not. 

The behavior of (t) has recently been obtained analyt- 
ically in the thermodynamic limit in the quench setting 1 17|. 
The results of iflTl for {t) can be summarized as follows. 
First of all (t) = identically for quenches starting in the 
disordered phase {\hi\ > 1) as the symmetry remains unbro- 
ken. For quenches starting in the ordered phase {\hi\ < 1) 
one has a different behavior depending whether one ends up 
in the ordered or disordered phase. More specifically 



(t) = 



Ae-*/^ I/12I < 1 

A'e-*/^^/! + cos {Ak„t + a) + ... \h2\ > 1 



The constants A, A' , r, /cq, ct all depends on hi, /i2 and are 
given explicitly in ifTTl . For quenches starting end ending in 
the ordered region, the equation (t) = = (00) = 
has no real solution. Correspondingly at finite size, T™ 
must be an increasing function of L and the simplest guess 
is T™ (X L'' — L. Instead for quenches ending in the dis- 
ordered region, we see that (<) = has a finite solution 
even in the thermodynamic limit, and so we expect in this 



case T' 



cq 



III. THE ANDERSON MODEL 

We now turn to random systems. The model we consider is 
the tight-binding model with random diagonal disorder, some- 
times referred to as the Anderson model lU : 



H 



L 

E 



(7) 



The diagonal elements /i^ are identically distributed indepen- 
dent random variables. In all of our simulations we will use 
a uniform, flat distribution in the interval [—W, W]. Through 
the Jordan- Wigner mapping, Hamiltonian Eq. (|7]i equivalently 
describes an XX chain of L spins in a random magnetic field. 
The Hamiltonian Eq. dTli can be written in compact notation 
as H = cJ Mf^iC , where c = (ci, . . . , cl) , Af^ is the one 
particle Hamiltonian and the subscript fi refers to the random 
variables /i^. In the infinite volume limit, the spectrum of il/^ 
is given, with probability one, by a (V^) + supp (— /^i) = 
[— 2i, 2t] + supp {—fJ-i) where is the discrete Laplacian 
in ID, (V^)^^. = Sij+i + Sij-i, and supp is the sup- 
port. In case of a uniform distribution we have a (A/^) = 
[—2t ~ W,2t + W]. Moreover, for any finite amount of ran- 
domness, the spectrum is almost surely pure point, i.e. con- 
sists only of eigenvalues, and the eigenfunctions are expo- 
nentially localized (see e.g. i fTSl IT9l ). Such a situation is re- 
ferred to as a localized phase. Since in the absence of disorder 
the model Eq. (j?]) is a band conductor, in one spatial dimen- 
sion there is a metal-insulator transition for any however small 
amount of disorder > 0. 

To study the equilibration properties of the Anderson 
model, we proceed as follows. We initialize the system in 



a state po, evolve it unitarily with one instance of Hamilto- 
nian Eq. (j?]) into p (t), and consider the expectation value of 
some observable A: A {t) — tr {Ap {t}). For very large sys- 
tems one expects that concentration results will apply and that 
the single instance A (t) will be, with very large probability, 
close to its ensemble average A (t) := E[^(/;)] where we 
denoted with E [•] the average over the random potentials /i^. 
When this happens, or more precisely when the relative vari- 
ance A^A/A^ as the system size increases, one says that 
A is self-averaging. We will not be concerned with this issue 
here, we just notice that A {t) generally gives the result of a 
hypothetical measurement of A within the confidence inter- 
val. The equilibration time is then obtained by solving for the 
first solution of A(Teq) = A. 

To be concrete we will consider a Gaussian initial state po 
with N particles uniquely specified by the covariance matrix 
= tr(c|cj/9o)- Since Eq. (g conserves particle number, 
the evolution is constrained to the sector with N particles. As 
observables, we choose a general quadratic operator given by 



In this case the expectation value A (t) 



can be completely characterized in terms of one-particle ma- 
trices |20|: 



A (t) = tr [Ap {t)] = tr (c 



(8) 



To be specific we will study two particular quadratic observ- 
ables: :— Yfj^i'^l^i which counts how many particles 
are present in the first £ sites (say from left). In this case the 
thermodynamic limit is given by fixing the particle density 
h' = N/L together with the "observable" density a — £/L 
and let L 00. 

As previously argued, a useful quantity to consider is the 
LE. In this free Fermionic setting it can be written as 11211 l22ll 



C{t) = |det(l-i?^ + i?^e 



)l 



Because of the random nature of the problem we do not 
expect the initial locations of the N particles to matter par- 
ticularly. Hence we will consider an initial state where all 
the N particles are pushed to the left, i.e. Ri i — 1 for 
i — 1,2, ... ,N and all other entries zero. We have checked 
that initializing the particles at other sites does not change 
our results. In any case for pure initial states, R-^ is a (or- 
thogonal) projector (i?^) — R^, meaning its eigenvalues 
are either 1 or 0. Therefore, in some basis, R^ always has 
the aforementioned form, i.e. there exists a unitary matrix X 
such that X^R^X = R% = diag (1, 1, . . . , 1, 0, . . . , 0) with 
N entries 1 and {L — N) zeros. The Loschmidt echo then 

becomes = \det {l - Rj^ + R%X^ er'*^'' X)\^ . Our 

choice of initial state corresponds to having X — 1. Thanks 
to the simple form of i?^, one can evaluate the determinant 
using Laplace's formula and reduce it to a determinant of an 
N X N matrix, i.e. 



£{t) = jdetr^v [X^e-"*^ 



where the operator Fat [Y] truncates the last L — N rows and 
columns of Y. Since clearly UFat [[/] ^|| < ||^|| for any uni- 
tary matrix U and L-dimensional complex vector ^, the eigen- 
values of Fat [X^^e^***^X] , Zi (t) have modulus smaller than 
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Figure 6: Relaxation of ne (t) — E [Ne (t) /L] with time for lengths 
L = 16, 40, 64, 112, at filling ly = N/L = 0.25 and a = £/L = 
0.5. The random distribution has variance a = 0.6 (W — 0.6\/3 ~ 
1.03). The ensemble average is computed using 1000 realizations at 
each time. 

one. The LE can than be written as 

C (t) = e^(*) = e''^/^*) 

having defined / (t) = (1/^^) E^i In k» Wl^- Provided 
/ (t) has a limit as L oo, this shows that C (t) is exponen- 
tially small in the system size. Accordingly, since J- {t) := 
In C (t) is extensive, we expect it to be self-averaging and con- 
venient to study. We also consider the averaged Loschimdt 
echo L{t) — E [C (t)] and the ensemble average of the loga- 
rithm of the Loschimidt echo (LLE) F (i) = E [T (t)]. Note 
that by Jensen's inequality L (t) > e^^^\ 

In Fig. |6] we show the results of our numerical simu- 
lations for the -ensemble averaged- observable (t) = 
E [Ni (t) /L] . When computing the equilibration time for ob- 
servable A by looking for the first solution of A(Toq) = A, 
the computationally most demanding part is the calculation 
of A, especially for large system sizes. For quadratic ob- 
servables, this computation can be simplified considerably. 
We first note that the time and ensemble averages com- 
mute. This is essentially a consequence of Fubini's theo- 
rem and the fact that all our quantities are bounded. The 
time average of A (t) for a particular realization can be com- 
puted diagonalizing M^. With the notation AI^ ~ Q^AQ, 
A = diag {ei, . . . , e^}, and Q unitary, one has A{t) — 

E„,™ (Q^^0)„,™ iQ'R^Q)rn,n <r^'^'--'-) . NoW, sim- 

ply note that, with probability one, the spectrum e„ is non- 
degenerate. This implies that the time average is given exactly 
by 

A = Y.{Q'AQ)^^^{Q'R'Q)^^^ (9) 

n 

The value A is then computed taking the ensemble average 
of ^, A = E \A\ , using Eq. Jol. In figurejrjwe plot the equili- 
bration time obtained for (T) as a function of system size L. 
Our numerical results show that the equilibration time scales 
exponentially in the system size. 
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Figure 7: Scaling of equilibration time Tcq{L) vs system size L, 
extracted from ne (t). Tcq is obtained from this first solution of 
E[nf(rcq)] — E[ne{t)]. Parameters are as in Figure [5] . The en- 
semble averages are obtained summing over 2000 realizations. The 
equilibration time is found to obey log(req) = 0.021L°-^^^'' + 2.53 
with residual norm 0.17. 

For the case of the LE and the LLE we have not been able to 
compute the time average exactly. It is known that in case of 
non-degenerate many body energies, which is safe to assume 
in presence of randomness, one has C — ^„ | {tpo \ {n}) |^ lH 
(where | {n}) — \ni, . . . , 7ii)'s are the many-body eigenfunc- 
tions satisfying J^j "-j — For our choice of initial state, 
the amplitudes are given by 

(V'olM) -dety[i^2,...,iv]Q^^/„} 

where V[i jv] is the N x L matrix with ones on the diagonal 
and zero otherwise. is formed in the same way but the 
ones on the diagonal are in correspondence of the row i for 
which Hi = 1. The normalization of the weights is provided 
by Cauchy-Binet's formula 

= detV[i^2,...,iv]Q^QV'[j2,...,^v] =1 
The time average Loschmidt echo is then given by 

|detl/[i,2,...,w]g^l/{^„jp 

The above sum, however, contains an exponential num- 
ber of terms and is not practical for numerical com- 
putations. To evaluate L and F, we compute the /x- 

random average l{t) iV-J^p^^, ^t"^ A (t) (and 
similarly for F) using as many as A^sampios = 2000 
for A'timcs random times uniformly distributed between 
[Tq,T]. The time average is then obtained via L ~ 

KlpiesN^iL. (T To)-' EjlT Etl"""^ ^^it,)■ We 
use a positive Tg to get rid of the initial transient and ob- 
tain more precise estimates with the same computational cost. 
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Figure 8: Relaxation of F{t) — E[log{£{t))]. Parameters are as in 
Figure[6]and an ensemble average was performed over 1000 realiza- 
tions at each instant in time. 
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Figure 9: Scaling of the relaxation time Tg, such that F{Teq) — F . 



Te, = 1.39L" 



0.847 with residual norm 0.05. 



Typical valued we used are T ~ 1000, Tq ~ 750. Larger 
times of the order of T ~ 10000 with Tq ~ 1000 were used 
for systems of size exceeding L ^ 128. 

In figure [8] we plot the ensemble averaged LLE time series 
F (t) at different sizes L while in figure |9] we show the equi- 
libration times of the logarithmic Loschmidt echo as a func- 
tion of L . In this case as well, the numerical results indicate 
equilibration times diverging exponentially in the system size. 
Based on the numerical evidence, we conjecture that for more 
general observables O, the equilibration time in the Anderson 
model might satisfy 

where cq, do, ip are constants which depend on the observ- 
able O. 



IV. A CUE MODEL 

In this section we consider another random matrix model 
for which we are able to prove analytically that Tcq < 1. In 
this model we fix the evolution operator U at time t = 1 
to be taken from the circular unitary ensemble (CUE). This 
means that U is an L x L unitary matrix sampled from 



F(t) (ti ± a) 
- Asymptotic 



l-U 




Figure 10: Plot of F {t) for L = 120, at = 1/2 (half tilling). The 
variance 5F (t) is shown as error bar. Numerical results are obtained 
averaging over A^sampics = 7500 unitaries distributed according to 
the Haar measure and obtained with the algorithm outlined in 1231 . 



the uniform Haar measure over the group U (L). At other 
times the evolution is defined via [/*. The arbitrariness in 
the definition of [/* for t e M is fixed in the following 
way. Any unitary matrix U from the CUE can be written 



as [/ = V^e^'^V where V is again Haar distributed, 
diag(e 



,i4> 



i(j}2 



g#i.) and the phases (pi are distributed 



ai't'i I 



where C ^ ~ 
e [0, 27r) (see 



according to P (0) C'ni<j 
(27r) L\ is the normalization constant and 
e.g. l|24|). Our model is described by taking a Hamiltonian 
H — 4> := diag((/)i, (f>2, ■ ■ ■ , <j>L) and considering the aver- 
age over all isospectral Hamiltonians H' = V^HV with V 
Haar distributed. The dynamical evolution is given by e**^ . 
Given these considerations this model is equivalent to the en- 
semble considered in |6 , 7| after averaging all the energies Ei 
with the CUE distribution P {4>). Let us now turn to the com- 



putation of F (t) 



Indetr 



N 



lit 



integer, the evolution operator is given by e'*^ — where 
J7" is unitary and Haar distributed, hence, at integer times we 
obtain 



F in) = Fo 



ln|detrAr([/)|' 



independent of n. We observed, as it is natural to expect, that 
the function F (t) has a limit as i oo (see Fig. 10 1. In this 



case the limit must coincide with the time average and with 
F(?i) for n integer, i.e. F (oo) — f{n) — F. This shows that 
the relaxation time in this random systems is bounded by one. 
In fact in principle one could have f (t) ~ f also for a time t 
smaller than one, however our numerics indicates that the first 
occurrence of f (t) — F is indeed at Tcq = 1 (see Fig. 10 1 
independent of the system size. We would like to mention at 
this point the results of Ref. [1 1 where a different sparse ran- 
dom ensemble has been constructed and an equilibration time 
growing as the system size has been reported. These findings 
show that random systems can give rise, in general, both to 
slow and fast equilibration processes and the correct equili- 
bration time-scale can only be obtained through an accurate 
investigation (although one expects faster equilibration to be 
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associated to less sparse Hamiltonians). 

It turns out that the Umiting value F (oo) can be obtained 
exactly. The distribution of eigenvalues of truncated matri- 
ces rN{U) when U is Haar distributed has been computed 
in [25J. The eigenvalues zi of rN{U) are complex num- 
bers in the unit disk \zi\ < 1. Calling = and defin- 
ing the probability distribution of the moduli Pl.n {'>') = 
r- rr^__ ' the CUE average at integers time 



N~'rli^u [Sir 
is then given by 



fo^N [ Pl,n (r) In (r^) dr. 
Ja 



(10) 



Zyczkowski and Sommers were able to compute the distribu- 
tion of the moduli PL,N{r) and obtained, with x = r'^ 



L-JV-l 



P(r) 



2r (1 - x) 

N {L-N-iy. ydx 



d \ 1- x^ 



l~x 



The probabiUty density (r) in the thermodynamic limit at 
fixed particle density N/L = v, L ^ oo, was also computed 
in ||251 and is given by 



PAr) 



1 



2r 



(11) 



for r < and zero otherwise. Plugging Eq. (Ill into 



Eq. ( [T0| one obtains, in the thermodynamic limit, the follow- 
ing particularly simple result 



F (cx)) = i [i/ In 1/ + (1 - t/) In (1 - i^) 



(12) 



Quite surprisingly Eq. ( [T2| is the negative Von Neumann 
entropy of the Gaussian state with covariance matrix R (n) = 
E [R (n)] obtained taking the ensemble average of the covari- 
ance matrix at integer times t ^ n. The reasoning is the fol- 
lowing. First we remind that the von Neumann entropy of 
a Fermionic Gaussian state pw with covariance matrix W is 
given by 

- S'vN (pw) = trpvK Inpw = 

tvW In W + tr (I - W) In {1 - W) . 

Now note that at integer times the evolution operator is L/" — 
U' with U' again unitary. At such times the average of the 
covariance is then R (n) — Ejj [J7i?[/^] and is proportional 
to the identity by Schur's lemma. The constant is fixed noting 
that trR (n) = tri? = N, which implies R [n] = vl. The 



claim then follows trivially taking traces of diagonal opera- 
tors of the form (7I) ln(7l). We do not know the average 
covariance at non-integer times, however if a limit t — > 00 
exist it must coincide with R(n). In that case we would 
have F(oo) = -S (/Or(oo)) [= -S (pr(„))]- Research is in 
progress to check weather this connection between the aver- 
age of the logarithmic Loschmidt echo and the von Newmann 
entropy has more general validity. 

Finally, we also verified that F (t) is indeed self averaging, 
i.e. the relative variance goes to zero as L increases. In par- 
ticular, since F {t) cx L, the variance of the rescaled variable 
F [t) /L goes to zero. 

V. CONCLUSIONS 

In this paper we considered equilibration in finite- 
dimensional isolated systems, and particularly concentrated 
on the time needed to reach equilibrium and its scaling be- 
havior with the system size. The standard definition of equi- 
libration time extracted from the exponential decay of some 
observable does not work in finite system because the dynam- 
ics is quasi-periodic and thus, no exponential decay can take 
place. Our definition of equilibration time Toq is precisely the 
time needed for an observable A (t) = ("^ol^ (t) \ipo) to reach 
its equilibrium value, and is given by the earliest solution of 
A (t) ~ A. We first examined clean systems. Considering 
the Loschmidt echo as a particular observable, we showed 
that in the general situation of a gapped or clustering initial 
state, the equilibration time is independent of system size, 
i.e. Toq = O On the other hand, for small quenches 

close to a critical point, one finds Tcq = O (i^) where C, is 
the dynamical exponent. We then turned to random systems 
and tackled the tight-binding model with diagonal impurities 
as an example. Considering different observables A, we found 
in all cases, that In Toq = c^T'^'-* + dA , where CA,dA, ipA are 
observable-dependent constants. The exponential divergence 
of equilibration time, however, seems a general feature of the 
localized phase in this system. Finally we introduced a novel 
random matrix model (similar to the one considered in |6, 7|) 
for which we were able to prove that Toq = O (1). This shows 
that, obviously, a higher degree of randomization can help the 
system reach equilibrium faster. 
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the ARO MURI grant W91 lNF-1 1-1-0268, and by the Na- 
tional Science Foundation under Grant No. NSF PHYll- 
25915. 
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